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I.  BACKGROUND 

The  effectiveness  calculations  for  Javelin  have  evolved  greatly  in  the  past  two  years.  The 
Systems  Simulation  and  Development  Directorate  (SS&DD)  has  been  tasked  to  calculate  Proba¬ 
bility  of  Kill  (Pk)  based  on  output  from  digital  and  real-time  flight  simulations.  The  U.  S.  Army 
Research  Laboratory  (ARL)  at  Aberdeen  Proving  Ground,  MD  is  producing  the  lethality  data  for 
these  analyses.  The  results  and  data  that  are  produced  by  the  simulations  at  SS&DD  will  be  for¬ 
warded  to  the  U.  S.  Army  Material  Systems  Analysis  Activity  (AMSAA). 

Software  has  been  developed  within  SS&DD  to  manipulate  official  Pk  cell  data  generated 
at  ARL  and  to  calculate  a  Pk  for  Monte  Carlo  sets  of  terminal  conditions  from  digital  and  real¬ 
time  flight  simulations.  1  .ie  software  also  allows  an  analysis  of  these  sets  through  the  use  of 
graphical  and  database  techniques.  All  software  was  developed  on  Silicon  Graphics  computers 
and  written  in  the  **C”  computer  language. 

Meetings  were  held  between  the  SS&DD,  the  Javelin  Project  Office,  ARL,  and  AMSAA  to 
determine  how  to  best  fulfill  each  agency’s  requirements  while  representing  the  system  in  a 
proper  manner.  This  report  will  discuss  the  lethality  effectiveness  computations  methodology 
and  distribution. 

II.  ANALYSIS 

The  SS&DD  effort  to  determine  system  effectiveness  initially  consisted  of  an  analysis  of 
the  terminal  conditions  from  the  digital  Six-Degrees-of-Freedom  (6-DOF)  simulation  and  the 
lethality  cell  maps  developed  by  ARL.  These  lethality  data  were  provided  for  certain  missile 
body  angles  (body  angle  theta).  Each  theta  consists  of  a  particular  trajectory  angle  (gamma)  and 
an  inertial  angle  of  attack  (alpha)  combination.  Figure  1  schematically  depicts  these  angles.  The 
same  variations  were  noted  in  the  yaw  plane  (Fig.  2).  Numerous  flight  simulation  runs  have 
shown  that  the  angles  may  vary  as  much  as  5  to  10  degrees  from  run  to  run.  This  implies  that  a 
mean  shotline  may  not  correctly  represent  the  best  (or  most  correct)  shotline  and  improper  credit 
would  be  attributed  to  the  system.  Figure  3  shows  a  representation  of  what  would  happen  when 
varying  angles  of  attack  are  applied  at  the  terminal  point. 

Through  analyses,  it  was  determined  that  the  existing  data  did  not  adequately  represent  the 
missile  terminal  state  vector,  with  its  variances  of  angular  orientation.  These  angular  definitions 
are  important  for  the  calculation  of  lethality  maps,  since  threat  vehicles  can  exhibit  large  discon¬ 
tinuities  from  one  cell  to  an  adjacent  cell.  This  is  highly  dependent  on  the  location  of  critical 
components  and  whether  the  warhead  penetrated  (or  hit)  that  component.  Varying  a  few  centi¬ 
meters  in  component  location  or  degrees  in  shotline  orientation  can  make  a  large  difference  in 
the  Pk  of  a  cell.  It  should  be  noted  that  the  initial  data  request  was  developed  by  the  Javelin 
prime  contractor  to  obtain  design  points.  These  data  were  not  meant  for  use  in  a  system  perfor¬ 
mance  assessment.  For  this  reason  SS&DD  requested  a  high  resolution,  angular  oriented  Pk  data 
set.  This  data  covers  a  two  sigma  variation  of  the  terminal  angles,  in  five  degree  increments. 

To  calculate  Pk  for  the  Integrated  Flight  Simulation  (IFS)  and  Hardware-in-the-Loop 
(HWIL),  terminal  conditions  are  received  from  the  digital  IFS  and  HWIL  Simulation.  Data  re¬ 
ceived  from  these  simulations  include  a  three-dimensional  hit  point  (where  0,  0,  0  is  the  center 
of  the  turret  on  the  ground  plane),  target  azimuth,  the  trajectory  angle  in  the  pitch  plane  (yp),  tra¬ 
jectory  angle  in  the  yaw  plane  (yy),  inertial  angle  of  attack  (a),  and  yaw  (P).  The  conditions  of 
the  scenario  run  (countermeasures,  background,  time  of  day,  season,  etc.)  are  also  inputs  into  the 
evaluation  tools. 
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The  three-dimensional  hit  point  is  projected  into  a  two-dimensional  (h,v)  hit  point  in  the 
plane  normal  to  the  mean  trajectory  or  velocity  vector  (yp,  yy).  The  dispersions  (oh,  ov)  are  cal¬ 
culated  in  this  plane.  If  the  data  shows  a  dependency  on  season/time  (as  with  short  ranges,  or 
non-specification  cases)  separate  dispersions  are  calculated  for  each  season/time  within  a  Monte 
Carlo  set.  The  Pk  for  each  run  of  the  Monte  Carlo  set  uses  an  average  Pk  from  the  nine  cells 
(three  by  three)  closest  to  each  (h,  v)  hit  point  location  on  the  cell  map  selected  from  the  nearest 
terminal  angles  (yp,  yy,  a,  P).  This  will  be  the  Pk  for  that  specific  run.  These  Pk  values  are  av¬ 
eraged  together  to  determine  the  Pk  for  the  Monte  Carlo  set.  This  method  is  known  as  the  Local 
Pk  Method.  Probability  of  hit  (Ph)  is  calculated  by  dividing  the  number  of  hits  by  the  sum  of  the 
hits  and  misses.  Some  misses  are  statistical  outliers  that  indicate  some  form  of  system  failure, 
such  as  a  loss  of  lock  on  the  track  gates.  Since  these  outliers  are  system  errors  and  not  accuracy 
errors,  logic  to  exclude  them  from  the  (h,  v)  standard  deviation  calculation  has  been  incorpo¬ 
rated.  Outliers  are  defined  as  those  shots  impacting  outside  of  an  eight  meter  radius  from  the 
origin  of  the  tank  coordinate  system.  Two  types  of  Pk  are  defined,  the  Pk  given  a  shot  (Pks)  and 
the  Pk  given  an  engagement  (Pke).  Pke  includes  the  effects  of  the  system  rejecting  an  initial 
lock  of  the  target,  note  that  this  is  a  different  definition  than  that  defined  in  the  Javelin  Specifica¬ 
tions  Document.  The  definitions  are: 

Pks  =  ZPk  for  each  hit/  (#  hits  +  #  misses) 

Pke  =  ZPk  for  each  hit/  (#  hits  +  #  misses  and  rejections). 

Output  from  the  analysis  of  each  Monte  Carlo  set  included  the  identification  of  mean  termi¬ 
nal  angles,  Center  of  Presented  Area  (CPA),  mean  hit  point,  the  number  of  hits,  misses,  rejec¬ 
tions,  aborts,  and  outliers,  Ph,  and  Pks  and  Pke  for  mobility,  firepower,  mobility/firepower,  and 
catastrophic  kills.  Since  an  abort  is  a  simulation  failure  and  has  no  relation  to  the  performance 
of  the  system,  those  runs  are  excluded  from  any  calculations.  There  are  three  encoded  words 
that  show  the  conditions  incorporated  into  the  Monte  Carlo  set  (background,  location,  counter¬ 
measure,  etc.).  A  grid  of  the  target  will  be  provided  displaying  the  dispersions,  CPA,  Center  of 
Hitpoint  Dispersion,  and  hit  points  identified  by  the  four  different  season/time  variables.  Figure 
4  shows  a  sample  output  of  the  Monte  Carlo  analysis.  Several  types  of  files  were  created  for  fur¬ 
ther  analyses,  each  were  utilized  to  investigate  the  cardioid  or  Combat  Vehicle  Combat  Perfor¬ 
mance  Operations  Assessment  (CV-CPOA)  weighting.  These  files  are  created  for  use  at  other 
agencies.  One  such  file  is  the  SUMLIN  file,  a  sample  is  shown  in  Figure  5.  It  will  include  all 
runs  for  a  particular  target  or  topical  study. 

Since  there  is  such  a  large  amount  of  cell-by-cell  data  for  each  target,  there  is  some  diffi¬ 
culty  in  other  agencies  utilizing  the  data.  For  example,  AMSAA  required  a  single  cell  map  file 
for  each  range  and  azimuth,  to  perform  their  analyses.  SS&DD  created  a  single  file  by  taking 
the  cell  map  files  for  a  specific  target,  and  applying  a  normal  distribution  of  the  terminal  angles 
(yp,  yy,  a,  |3)  about  each  cell  on  a  plane  normal  to  the  mean  velocity  vector  at  that  range.  The 
normal  distribution  is  based  on  the  6-DOF  runs  completed  by  the  Joint  Venture  to  characterize 
the  systems  behavior.  The  equation  for  the  formulation  of  these  cell  maps  is 

Pk  =  Dyp  2a  Zyy  2p  Poccurrenceypayyp  Pkypayyp 

Poccurrenceypayyp  =  (fne(yp2/1.414)-fne(Ypl/1.414))* 

(fne(a2/1.414)  -  fne(al/1.414)  (fne(yy2/1.414)- 

fne(Yyl/1.414))  *  (fne(p2/1.414)  -fne(P  1/1.41 4))/l 6.0 
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where 


yp2  is  the  upper  bound  of  the  velocity  vector  in  the  pitch  plane 

Ypl  is  the  lower  bound  of  the  velocity  vector  in  the  pitch  plane 

a2  is  the  upper  bound  of  the  inertial  angle  of  attack 

al  is  the  lower  bound  of  the  inertial  angle  of  attack 

yy2  is  the  upper  bound  of  the  velocity  vector  in  the  yaw  plane 

yyl  is  the  lower  bound  of  the  velocity  vector  in  the  yaw  plane 

FNE  is  a  standard  error  function  found  in  most  statistical  manuals. 

Figure  6  shows  an  example  of  the  a  cell-by-cell  map  created  by  this  means.  Note  the  excess 
cells  on  the  edges  of  the  target,  these  contribute  very  little  Pk  and  may  be  ignored  when  applying 
an  analyses  to  the  CPA.  These  cells  are  caused  by  the  variations  in  the  occurrence  of  the  termi¬ 
nal  angles. 

in.  CONCLUSIONS  AND  RESULTS 

The  SS&DD  will  provide  data  for  the  overall  Javelin  Performance  Assessment  Any 
topical  studies  needed  for  system  evaluation  were  also  evaluated  with  this  effectiveness  method¬ 
ology.  Pk  and  accuracy  data  will  be  provided  from  the  IFS  and  HWIL  simulations.  Cell-by¬ 
cell  data  for  specific  ranges  that  are  independent  of  the  natural  angular  variations  will  be  created 
and  supplied  to  various  government  and  contractor  organizations.  This  methodology  will  be 
employed  for  all  targets  run  in  the  Javelin  Performance  Assessment.  Results  of  the  Performance 
Assessment  will  be  published  in  a  separate  report. 
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Theta  =  Gamma  Pitch  +  Alpha 


Gamma  Pitch 


(Front  View) 


Figure  1.  Missile  Body  Angle  Definitions  for  Pitch  Plane 


Psi  =  Gamma  Yaw  +  Beta 


Gamma  Yaw 


Azlmath 


(Top  View) 


Figure  2.  Missile  Body  Angle  Definitions  for  Yaw  Plane 
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Figure  4.  Sample  Output  for  Monte  Carlo  Analysis 
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Sample 


Kill  Type:  Mobility  /  Firepower 

Figure  6.  Cell-by-Cell  Map  Based  on  Weighting  of  Terminal  Angles 
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